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Abstract 



Nuclear reactions in stars occur between nuclei in the high-energy tail of 
the energy distribution and are sensitive to possible deviations from the stan- 
dard equilibrium thermal-energy distribution. We are able to derive strong 
constraints on such deviations by using the detailed helioseismic information 
of the solar structure. If a small deviation is parameterized with a factor 
exp{-5(E/kT) 2 }, we find that 5 should lie between -0.005 and +0.002. How- 
ever, even values of 5 as small as 0.003 would still give important effects on 
the neutrino fluxes. 
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A. Introduction 



Thermal averages are fundamental ingredients of the theoretical description of many 
physical phenomena: solar modeling is a specific example. These thermal averages can 
and are often described as integrals weighted by the appropriate equilibrium distribution 
functions; e.g., the mean square averaged velocity (v 2 ) of a particle in a gas is obtained from 
the integral of v 2 times the one-body velocity distribution function f(v), the pressure from 
the integral of pv times f(v,p) (p is the momentum), and so on [PJEJ. 

In the limit of non-interacting states, infinite volume and zero density, a single scale (the 
temperature or the average one-body energy) characterizes all the equilibrium distributions, 
which are described by the Maxwell-Boltzmann distribution (MBD). It is well-known that, 
even for non-interacting states, when the system is finite and/or the density is not zero, the 
distribution deviates form the MBD, and the resulting distribution (microcanonical, Fermi- 
Dirac or Bose-Einstein) is characterized by additional scales (total energy, Fermi energy, 
etc.). Similarly, the interaction could produce additional dynamical scales that modify the 
free distribution, e.g., 4 He is phenomenologically better described as a weakly-interacting 
Bose system than in terms of its fermionic constituents. In principle the thermal distribution 
of the effective weakly interacting degrees of freedom (the bosonic 4 He nuclei, in this exam- 
ple) could be dynamically calculated from the original strong interacting elementary parti- 
cles (the nucleons). However, theoretical calculations of thermal distribution functions for 
strongly interacting systems are very difficult, and one often resorts to phenomenologically 
motivated parameterizations. In specific cases, it has been possible to derive equilibrium dis- 
tributions departing from the MBD and, more in general, theoretical frameworks |^-|7j 
have been formulated that naturally produce nonstandard distributions. 

In spite of this, one can argue that, even in presence of strong and/or many-body and/or 
long-range forces, one single scale dominates in many practical case and, therefore, that 
the MBD is an excellent approximation. This argument is confirmed a posteriori by its 
phenomenological success. However, one should also keep in mind that many applications 
do not test the details of the distributions, but only one or a few moments. In particular, 
if a physical quantity is determined by only one moment, one can always summarize the 
relevant information in the most economical way by using the MBD. In practice, in many 
cases nothing changes using distributions that differ only in the higher moments. This low 
sensitivity of many important physical observables to the details of the thermal distribu- 
tion together with the difficulties of the microscopic calculation leads to the possibility of 
considering more general distributions that depart from the MBD. 

For instance, already two decades ago PHTOf it was proposed that small depletions of the 
high-energy tail of the relative energy distribution could modify the solar neutrino fluxes. 
This same idea has been recently reconsidered ||TTHT4|| in the light of the new developments 
that put nonstandard equilibrium distributions on a firmer ground. 

Similarly, one could invoke a small enhancement of the high-energy tail of the proton 
distribution in order to efficiently burn lithium near the bottom of the convective zone. This 
could be regarded as an attempt to account for the low photospheric lithium abundance 
(about a factor 100 lower than the meteoric one essentially unexplained within the 

Standard Solar Model (SSM) (see, however, Ref. [fTS] ). 

In this paper we take the opposite approach and study what kind of constraints our 
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best knowledge of solar physics, both theoretical and observational, can impose on possible 
deviations from standard thermal distribution. Our two basic tools will be the sub-barrier 
fusion reaction rates and helioseismology. 

Nuclear reactions in stars occur generally between nuclei with kinetic energy much larger 
than kT and are thus suitable for probing the high-energy tail of the particle distribution. 
Even for the pp reaction, which has the lowest barrier, the Gamow peak in the solar core is 
at energy about five times larger than kT, making the reaction rate very sensitive to changes 
in number of high-energy particles. Therefore, if we can precisely determine a reaction rate 
by means of observations, this determination can be used to constrain the particle energy 
distribution. Stellar interiors are indeed an ideal laboratory for this investigation: they 
are to a very high degree in thermal equilibrium, and the density is high enough to make 
deviations from standard statistics conceivable. 

Moreover, helioseismology allows us to look deeply in the core of the Sun. The extremely 
precise measurements of a tremendous number of frequencies enable us to extract values 
of sound speed with high accuracy even near the solar center [|17|--|22|| . In addition several 
properties of the convective envelope are accurately determined by means of helioseismology. 

Recent standard solar models that include the state-of-the-art "standard" solar physics 
are in good agreement with helioseismic data [|T^,^1|,|T^] . These solar models implicitly as- 



sume that the solar core can be described in terms of a gas of particles interacting via 
two-body forces with no many-body effects apart for mean-field screening. In particular, 
the ion relative velocity distribution follows the MBD and the diffusion of the average num- 
ber of particles is Brownian. In some sense, helioseismology tells us that this framework 
is basically correct. Nevertheless, it is important to quantitatively assess to what extent 
nonstandard distributions are still compatible with present data. In this respect we remark 
that the information on the solar interior provided by helioseismology is so detailed that the 
pp reaction rate can be reliably constrained ||23|| . 

We shall investigate solar models obtained by modifying a SSM so as to include the 
effects of a nonstandard distribution. By requiring that the predictions of the resulting 
solar models agree with the helioseismic determinations of convective envelope properties, 
we shall constrain the possible nonstandard distributions. 

B. Modified statistics and burning rates 

In the ordinary treatment, the single particle energy distribution is taken as a MBD for 
protons and other ions Q: 

fMBo(E) = -^= e~ E/kT ■ (1) 



The nuclear burning rate between nuclei with mass numbers i and j is given by: 



In the solar plasma quantum corrections to the statistics are relevant for electrons but not for 
nuclei. 
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(<™)ij = J dea ij (e) fij(e)v(e) , (2) 

where e is the collision energy, v is the relative velocity, the cross section a has the usual 
parameterization 

a«(6) = ^ e -^V^, (3) 

and fij(e) is the collision-energy distribution of the reacting nuclei. As well known, if the 
one particle distribution is a MBD, so it is the collision-energy distribution fij(e). 

Small deviations of /y(e) from the MBD can be parameterized to first approximation by 
introducing a dimensionless parameter S 

f^\e) = f MBD (e)e- S ^\ (4) 



so that for 5=0 the classical statistics is recovered |5,IU,rj|; for small 5 the distribution 
is close to the standard one at values of e near the thermal energy kT, whereas significant 
distortion occurs in the high-energy tail. For 5 > this parameterization implies a depletion 
of the tail. The same parameterization can also be used to mimic an enhanced tail (5 < 0), 



understanding that a suitable cutoff is introduced |L4 . 

Solar models corresponding to modified statistics have been built by using our stellar 
evolutionary code FRANEC ||24j| , where all the nuclear- reaction rates have been calculated 



according to Eqs. (Q) and We remark that we are assuming here that S is the same for 
every reaction and is constant in the nuclear energy production region (R/R & < 0.2). 
The solar structure is primarily sensitive to the rate of just two reactions: 

i) p + p — > d + e + + v e . Since this reaction is at the basis of the nuclear-reaction chain that 
sustains the Sun against gravitational collapse, it is natural that the internal solar structure 
is strongly influenced by its rate. As shown in Ref. f23fl , this rate is strongly constrained by 
helioseismic determinations of the convective envelope. 

ii) p+ 1A N — > 15 O+7. The rate of this reaction governs the efficiency of the CNO cycle, which 
is marginal according to the SSM. An enhancement of the high energy tail (5 < 0) makes the 
CNO cycle more efficient and even dominant, resulting in solar models drastically different 
from the SSM (e.g., the energy production is concentrated near the center, a convective core 
can arise, . . . ). 



C. Modified statistics and the properties of the convective envelope 



As shown in Ref. 



]17| , helioseismology determines with high accuracy three independent 
properties Q of the convective envelope: its depth R^, the density at its bottom p& and the 
photospheric helium abundance Y p h- 

The helioseismic values of these quantities are shown in Table | together with two esti- 
mates of their uncertainties: (AQ/Q) cons corresponds to the (very) conservative definition 
of Ref. |I7|], whereas (AQ/Q)i a is the corresponding la "statistical" error estimate. 

In the same Table we also show the predictions of the "model with helium and heavy ele- 
ments diffusion" of Ref. |25|] (BP95), which are in excellent agreement with the helioseismic 



4 



determinations. For this reason we shall use this model as the reference SSM. As an example 
of possible "systematic" theoretical uncertainties, we also show results from our solar model 
including helium and heavy elements diffusion [Ell (FR97), which deviates somewhat from 
the helioseismic determinations. 

By numerical experiments with FR97, we have also determined the dependence of these 
three properties on 5; results are shown in Fig. |I[ The different behavior for negative and 
positive 5's becomes more evident as \5\ increases. This qualitative difference reflects a 
physical effect: when the tail of the distribution is enhanced (5 < 0) the CNO cycle becomes 
important (at 5 = —0.01 the contributions of the pp-chain and of the CNO cycle are about 
the same). Since the Gamow energy for the p + 14 N reaction near the solar center is about 
27 keV, a factor five larger than the one for the pp reaction, even a small 5 yields drastic 
effects on the reaction rate and, consequently, on the solar structure. 

Nonetheless, for small values of 5, these dependences can be parameterized by power 
laws: 

tJSm = ( e Ql ' ( 5 ) 

where the constant exponents aQ i are shown in the last column of Table |. The solid curves 
in Fig. [I] demonstrate the goodness of such a parameterization in the range of 5 that is 
relevant to our results (\6\ < 0.005). 

Our basic strategy will be the following: we determine the acceptable range of 5 such that 
Rb, pb and Y p h are predicted within their helioseismic ranges, by using Eqs. ([5]) to determine 
the dependence of these properties on 5. 

There are at least four major uncertainties in building standard solar models that also 
have the potentiality of affecting the three helioseismologic properties under investigation, 
and, therefore, that could interfere with/hinder the effect of 5: the astrophysical factor S pp , 
the solar opacity k, the heavy element abundance ( = Z/X, and the diffusion coefficients. 
We shall add all these effects one after the other, and determine a range of #'s that takes 
into account these uncertainties. 



D. Results 

For determining the range of 6 allowed by helioseismology, we use several approaches 
corresponding to increasing conservativeness. We start by defining a \ 2 as: 

where Qi(5) are computed by using Eqs. @ and the errors are the la estimate of Table |. 
The value x 2 (0) indicates how well the SSM reproduces these helioseismic properties. The 
first row of Table |I| shows the good agreement between BP95 and helioseismology (x 2 /dof 
= 8.61/3). 

If we use 6 as free parameter (second row of Table H), we find the following best fit value 
( x 2 /dof = 0.08/2) and la range: 
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5 = (-0.77 ±0.26) x 10~ 3 



(7) 



These strict constraints on the allowed values of S come mainly from the precise determina- 
tion of density at the bottom of the convective envelope combined with its strong dependence 
on 5. In fact, the relative change of Qi due to A5 is approximately A5 x agr, therefore, the 
allowed variations of S can be estimated as A5 ~ (AQi/ Qi)x a /\aQ t \ (last column of Table |). 



1. Uncertainties on S pp 

A conservative estimate of the uncertainty is provided by the range of the published 
results pUfl , whereas a la estimate has been provided in [^7[; we shall use AS PP / S pp = 
0.05/3 at la (5% is the "3cx error" estimate). The dependence of Qi on Spp has been 
determined numerically in Ref. pfl . By redefining a suitable x 2 (A S pp ): 

X ^ Spp) ~ V { ) { *S PP ) ' (8) 

we find that the best fit value of S does not change and that the la range is double: 

5 = (-0.77 ±0.50) x 10" 3 , (9) 

and that, consistently, the best fit value for S pp is S^ssm- The facts that the SSM is already 
in very good agreement with helioseismology and that the dependence of the Qi on S pp is 
much weaker than that on 5 explain these results. 



2. Uncertainties on k and £ 



The heavy element abundance ( and the solar opacity k are known with a conservative 
accuracy of about 10% |28|j25|j29|] . Therefore, our la relative error estimate will be 0.1/3. 
The dependence?] of Qi on n and ( has been determined numerically in Ref. [23|. In this 
case, the relevant y 2 is: 



X 2 (5,S PP} (,k) = 



Qi(b~i Spp, £, k) Qd 



Qi 



Q — Q 



AQ t 
\ 2 

pp,SSM \ 



AS, 



K ~ KSSM_ V 

Ak ) 



C — CsSM \ 



(10) 

(ii) 



We find a small change of the best fit value and, again, an increase of the la range of 5: 

5 = (-0.75 ±0.67) x 10~ 3 . (12) 

Comparing the 3th 4th and 5th row of Table ||, one can notice that most of the effect is due 
to C- 



2 We remark that we are considering a constant rescaling of opacity along the solar profile. 
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3. Uncertainties on diffusion coefficients 



We use a SSM that includes element diffusion, calculated by solving the Burgers equa- 
tion |JtJ. Indeed, diffusion has been an essential ingredient of stellar evolutionary codes for 



achieving agreement between predicted and helioseismic values of properties of the convec- 
tive envelope |17| . The success of solar models with diffusion, and the corresponding failures 
of models that neglect diffusion, suggest that the diffusion process has been properly treated. 
However, in spite of the extensive discussion about the many assumptions underlying the 
calculation method [|31], j30| , p5|l , no quantitative estimate of the uncertainties of the calculated 



diffusion coefficients has been presented. 

Therefore, we also allow the diffusion efficiency to vary freely by rescaling the diffusion 
coefficients by an overall constant factor D (D = 1 corresponds to the SSM). We have 
determined the appropriate scaling laws [^2| of the properties of the convective zone. For 
completeness, we report the complete dependence on all the considered quantities (5, S pp , 
C, and D): 



Rh 



Rb,SSM 

Pb 
Pb,SSM 



Y, 



ph.SSM 




0.046 

0.016 



IT - 016 (13a) 
(13b) 
(13c) 



No additional term is added to x 2 , since we assume that D is completely undetermined 
(infinite error), and we let it vary freely. The only dependence of x 2 on D is through 
Qi($, S pp , C, n, D). As it is shown in the last row of Table ||, the la allowed range becomes: 

5 = (-0.91 ± 1.06) x 10" 3 , (14) 

and we find that the best fit value for D is only 3% smaller than the standard one. 



4- Solar model "theoretical uncertainties" 

At last we try to estimate how much our results could depend on having used BP95 as 
reference standard model. To this end, we consider one of the standard solar models (models 
that include all the state-of-the-art solar physics), whose helioseismic properties differ the 
most from BP95 and, consequently, fit less well the experimental data. We repeated the 
above- described analysis by using FR97 as standard solar model. When all parameters are 
varied, the la range and best fit value, cf. Eq. fli~4]), become: 

5 = (-1.79 ± 1.04) x 10~ 3 . (15) 

The corresponding fit to the helioseismic properties is acceptable (x 2 /dof = 2.32/1), the 
best fit value for D is 9% smaller than the standard one, and the values for ( and k are 2% 
larger. 
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E. Discussion and conclusions 

The constraints we have found on non-Maxwellian statistics look rather strict, as the di- 
mensionless parameter cannot exceed a few per thousand. In fact, if we define a conservative 
interval as the union of the 3a ranges found by using BP95 and FR97 SSMs, we find 

-4.9 x 10" 3 < 5 < 2.3 x 10" 3 . (16) 

However, even these small values of 5 could have non-negligible implications for those 
observables that are sensible to the high-energy tail of the distribution. As an example, 
we have estimated the possible effects of 5 = ±3 x 10 -3 in the two cases mentioned in the 
introduction. 



1. Neutrino fluxes 

In Table [II], we report the effect of nonstandard statistics on the main fluxes and on the 
signals of the chlorine and gallium radiochemical experiments. Even for such small values 
of 5 the boron and beryllium fluxes change substantially. 



2. Lithium abundance 

As well known, the photospheric abundance of lithium is a factor about 100 lower com- 
pared to the meteoric one fll5| . Different mechanisms have been proposed to explain this 



depletion [|T6| , p3| , |34| . Let us discuss the possibility that nonstandard velocity distribution 
could contribute to this depletion. First of all we note that, since the lithium abundance 
should be reduced in order to solve/alleviate the problem, the lithium burning rate should 
be enhanced relative to the standard case. This is achieved by a longer high-energy tail, 
i.e., 5 < 0. 

We assume that the limits on 5 derived in the production region apply also up to the 
bottom of the convective zone, and consider 5 = —3 x 10~ 3 . This value of S yields a reduction 
of the 7 Li abundance by only 7%, where the characteristics of the bottom of the convective 
zone has been taken from FR97 (T b = 2.1 x 10 6 K, p b = 0.18 g/cm 3 and X = 0.744). 
Depletions comparable with the observed ones could be obtained with 5 ~ —0.15, a value 
well outside the range reported in Eq. (|i~6|). 
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TABLES 



TABLE I. The three independent properties of the convective envelope used in our analysis. 
The first column (Qi) labels the property, the second (BP95) and third (FR97) columns show 
the values predicted by the reference solar model BP95 |^5[ and by FR97 |24| standard model, 
the fourth column (Qiq) shows the value derived by helioseismic measurements, and the next two 
columns the corresponding conservative and la errors. The last two columns show the exponents 
that determine the dependence from 5, aQ i = —dlogQi/dS, and the ratios between the la error 
and IckqJ (see text). 





BP95 


FR97 


QiQ 


(AQj/ Qi)cons 


(AQi/Qi) la 


a Qi 


(AQi/QOia/laqJ x 10 3 


Yph 


0.24695 


0.2321 


0.249 


0.042 


0.014 


6.2 


2.2 


Rb/R& 


0.712 


0.715 


0.711 


0.004 


0.002 


-2.2 


0.9 


p b [g/cm 3 ] 


0.187 


0.182 


0.192 


0.037 


0.0094 


33.6 


0.3 



TABLE II. Deviations from standard statistics allowed by helioseismic measurements. The 
first five columns show whether the parameter is kept fixed (F) at its SSM value or it is allowed to 
vary (V) as a free parameter within the range discussed in the text. The sixth column shows the 
resulting x 2 P er degree of freedom. The last two columns show the best fit value for 5 and its la 
error. 



D K 


C 


Spp 


5 


X 2 /dof Si 


Best X 10 3 


AS 


F F 


F 


F 


F 


8.61 /3 






F F 


F 


F 


V 


0.08 /2 


-0.77 


0.26 


F F 


F 


V 


V 


0.08 /2 


-0.77 


0.50 


F F 


V 


V 


V 


0.08 /2 


-0.73 


0.67 


F V 


V 


V 


V 


0.04 /2 


-0.75 


0.67 


V V 


V 


V 


V 


0.001/1 


-0.91 


1.06 


TABLE III. 


Effects 


of nonstandard statistics on neutrino fluxes. Relative deviations 


from SSMs 


of the Be and 


8 B neutrino fluxes and of the expected 


signals for gallium and chlorine detectors in 


two nonstandard solar 


models with deformed velocity distribution (<5 ^ 0). 












S = +3 x 10~ 3 




5 = 


-3 x 10~ 3 


A$Bc/$B e 






-0.30 






+0.38 


A$ b /$b 






-0.55 






+1.15 


A5ci/<Sci 






-0.16 






+0.31 


AScia/Sca 






-0.50 






+1.03 



11 



FIGURES 



FIG. 1. Dependence on 5 of the three independent properties of the convective envelope used 
in our analysis. Crosses show the values of the photospheric helium abundance Y p h (a), the density 
at the bottom of the convective envelope pb (b), and the depth of this envelope Rb (c), relative to 
their standard values, as functions of 5. The solid curves are the fits in Eq. (||) with the exponents 
otQi from the second last column of Table |. 
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